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Scope of NASA-Ames Cooperative Agreement No. NCC 2-361 


Work under this agreement started on August 1, 1985. This final report is pre- 
sented in two parts. The first part refers to correlation between prediction 
based on a linear theory and the U.S. Army's 7X10 ft. wind-tunnel data on flap- 
lag stability in hover and forward flight. The second part refers to the 
effects of quasisteady stall effects on the correlation. In accordance with the 
Cooperative Agreement Special Provision, the University and Ames personnel 
collaborated on performance of research, and the papers resulting from the 
research are co-authored. The papers are: 


1) Gaonkar, G.H., McNulty, M.J. and Nagabhushanam, J., "An Experimental and 
Analytical Investigation of Isolated Rotor Flap-Lag Stability in Forward 
Flight", 11th European Rotorcraft Forum, London, England, September 10-13, 1985, 
Paper No. 66. 

2) Nagabhushanam, J., Gaonkar, G.H., and McNulty, M.J., "An Experimental and 
Analytical Investigation of Stall Effects on Flag-Lag Stability in Forward 
Flight", 13th European Rotorcraft Forum, Arles, France, September 8-11, 1987, 
Paper No. 6-2-2. 


Part 1 - Linear Theory 


ABSTRACT 


For flap-lag stability of isolated rotors, experimental and analytical 
investigations are conducted in hover and forward flight on the adequacy of a 
linear quasisteady aerodynamics theory with dynamic inflow. Forward flight 
effects on lag regressing mode are emphasized. Accordingly, a soft inplane 
hingeless rotor with three blades is tested at advance ratios as high as 0.55 and 
at shaft angles as high as 20°. The 1.62-m model rotor is untrimmed with an 
essentially unrestricted tilt of the tip path plane. In combination with lag 
natural frequencies, collective pitch settings and flap-lag coupling parameters, 
the data base comprises nearly 1200 test points (damping and frequency) in for- 
ward flight and 200 test points in hover. By computerized symbolic manipula- 
tion, a linear analytical model is developed in substall to predict stability 
margins with mode indentification. To help explain the correlation between 
theory and data it also predicts substall and stall regions of the rotor disk 
from equilibrium values. The correlation shows both the strengths and 
weaknesses of the theory in substall ( | angle of attack | s 12° ). 
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azimuth position of the first blade). 
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Rotor shaft angle, positive nose down 
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Equilibrium pitch angle= 8 o + 0 s sin^j-+8 c cost!} c 
Total inflow ratio * ptana s +v 

Steady induced flow = v Q +v s r sinm^ + v c r cos^ 

Uniform, side-to-side and fore-to-aft unsteady 
(steady) induced flow 

Rotor rotational speed in rpm 

Dimensionless (1/ft) uncoupled lag frequency 

Azimuth position of the k-th blade 

Lock number 

d/dt 

absolute value 
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1. INTRODUCTION 


Relatively few investigations with concomitant corroboration of test data 
have been conducted on isolated rotor instabilities in forward flight. 

Such an investigation is particularly desirable for lag modes which are at best 
weakly damped and for which the state-of-the art of predicting damping levels 
merits considerable refinements.^ - ^ Here, we study a basic aspect of that 

investigation the adequacy of a widely used quasisteadv aerodynamics theory 

in predicting the flap-lag stability in forward flight (0_<p<^.55). 

We use an experimental model which closelly approaches the simple 
theoretical concept of a hingeless rotor as a set of a rigid articulated 
flap-lag blades with spring restraint. This intentionally built-in structural 
simplicity facilitates isolate the aerodynamic ingredients that participate in 
the correlation. The analytical model is based on quasisteady aerodynamics with 
dynamic inflow and it is directly tailored to fit the experimental model. The 
correlation refers to the lag regressing mode ^ which, compared to other modes, 
is the lowest-f requencv mode and is better suited to assessing the validity of 
the quasisteady aerodynamics theory. Although the theory merits substantial re- 
finements in stall ( |«| > 12'’' ), we have included some correlation in stall as 
well. Such inclusion provides an improved interpretation of the correlation 
between measured and predicted data. 
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A specific three-bladed rotor is investigated under untrim conditions^ * 1 1 
( e s = 0 c =O, cyclic flapping present) and it is soft in lag for forward flight 
test cases. However, from an intentionally simplified experimental model, the 
correlation makes use of a comprehensive data-base of nearly 1400 test points 
with respect to Q, R, 9 0 , a s and p and includes cases that vary from virtual 
zero-thrust conditions to aerodynamically demanding cases of highly stalled 
conditions with advance ratio p as high as 0.55 and shaft-tilt angle a s as high 
as 20". Such a correlation should give generally applicable qualitative results 
on the adequacy of the linear theory when stall is not an issue and should 
promote further refinements. 


2 . EXPERIMENTAL MODEL 

The model tested was a three-bladed hingeless rotor with a diameter of 1.62 
m. The rotor was designed as a set of rigid, articulated blades with spring 
restraint and coincident flap and lag hinges. This was accomplished by the 
"folded back" flexure design shown in an exploded view in Fig. 1. Taking the 
center of the flexural elements as the effective hinge point gives a 
non-dimensional hinge offset. of 0.11 for the design. The blades were designed 
to be very stiff relative to the flexures so that the first flap and lag modes 
involve only rigid body blade motion. The influence of the torsional degree of 
freedom was minimized by keeping the first torsion frequency as high as 
possible. In particular, the first torsion mode frequency was above 150 Hz 
non-rotating, insuring a rotating first torsion frequency of at least 9/rev over 
the entire rotor speed range tested. The rotor properties are summarized in 
Table 1. 
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The rotor had no cyclic pitch control, and collective pitch was set 
manually. The blade pitch could be set by changing the angle of the blade 
relative to the flexure at the b lade- flexure attachment, giving a structural 
coupling value of zero, or by changing the angle of the entire blade-flexure 
assembly relative to the hub, giving a structural coupling value of one. 

The test stand on which the rotor was mounted included a roll gimbal which 
could be locked out mechanically. Rotor excitation was accomplished through 
this gimbal by a 50 lb (222.4N) electromagnetic shaker. The stand was designed 
to be as stiff as possible so that the test would closely represent the case of 
an isolated rotor. However, the lowest frequency of the installed stand with 
the gimbal locked and the rotor mounted was found to be 31 Hz, somewhat lower 
than what was desired. The entire test stand could be pitched forward with an 
electric actuator to control the rotor shaft angle. The shaft angle provided 
the only means of controlling the rotor loads at a given collective pitch and 
advance ratio. A photograph of the model installed in the wind tunnel is shown 
in Fig. 2. 


3. ANALYTICAL MODEL 

The analytical model consists of an offset-hinged rigid lag-flap model with 
flap and lag spring restraints at the offset hinge. The spring stiffnesses of 
the math model are selected such that the uncoupled rotating flap and lag 
natural frequencies coincide with the corresponding first-mode non-rotating 
natural frequencies of the experimental rotor blade. The effect of the hinge 
offset, which is 11.1a in lag and flap, is accounted for in rotor trim and 
stability analyses. The rotor is untrimmed, with an essentially unrestricted 
tilt of the tip path plane (9 S = 8 C *0, and cyclic flapping present).^’ H The 
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rotor angular velocity is Q, and with time unit l/£2, the azimuth angle of the 
first or the reference blade represents the dimensionless time t. The 
distribution of flexibility between the hub including flexures and the blade 
(outboard of the blade location where pitch change takes place) is simulated by 
an elastic coupling parameter R which is also referred to as the hub rigidity 
parameter. Basically, R relates the rotation of the principal axes of the 
blade-hub system and the blade pitch 0. 

The blade airfoil aerodynamics is based on linear, quasisteady theory in 
substall (angle of attack | ct | <12 1 ') without inclusion of compressibility or 

other effects due to reversed and radial flow. Steady non-uniform inflow with 
first - harmonic distribution ( v ■ v 0 + v s r cosily + v c r sin liq^is assumed. The 
airfoil nonlinear section effects are neglected . ^ ^ For the math model we take 
a=5.73 and cj * 0.0079 ( NACA 23012). ^ ® The dynamic inflow effects are 

included from a 3x3 inflow model based on an unsteady actuator-disk theory the 

Pitt Model, which has been recently verified experimentally. ^ The equations of 
motion including the multiblade coordinate transformation are derived from 
symbolic manipulation 12,13. For the three-bladed rotor, the 15x1 state vector 

comp^rises 12 multiblade components (B 0 > £ 0 , B s , B* s , B c , B* c , Co* C* 0 » <f s < ?c> ?*c ] 
and 3 dynamic inflow components (v Q , v s , v c ). The computational details of 
Floquet transition matrix, and of damping levels and Floquet eigenvectors are 
respectively as in references 14 and 13. 
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4 . EXPERIMENTAL PROCEDURES 


An experimental investigation of the rotor's lag stability characteristics 
in hover and, especially, forward flight was conducted in the Aeroflight 
Dynamics Directorate's 7-by-10 foot (2.1-by-3.0 m.) wind tunnel. The data was 
collected in two tunnel entries, the first in the summer of 1982 and the second 
in the summer of 1983. The model's configuration and the test procedures used 
were the same for both tests. The only exception to this is that for the 1982 
test the hover data was taken with the wind-tunnel test section doors open and 
the windows removed in an attempt to reduce recirculation effects. In 1983 the 
hover data were taken incidental to forward flight data (ft=750 and 1000 rpms ) so 
the windows were left in place and doors were closed. The rotor plane itself was 
located about 0.63 rotor diameter above the test section floor (with the shaft 
vertical). However the influence of recirculation and ground effect on the 
hover data are not significant. ^ 

For each data point obtained, the blade pitch was set manually and the 
rotor was tracked. The rotor was then brought up to speed and wind tunnel 
dynamic pressure was increased to obtain the desired advance ratio while 
adjusting the shaft angle to keep the rotor flapping loads within limits. Once 
at the test condition, the roll gimbal was unlocked and the shaker was used to 
excite the model at the appropriate frequency. When a good level of excitation 
was evident in the output of the lag bending gages, the shaker was cutoff, the 
roll gimbal locked out, and then transient data acquisition commenced. The data 
itself was digitized on-line with a sample rate of 100 Hz and a total record 
length of 5.12 seconds. The data was then transformed to the fixed system 
coordinates. Spectral analysis and the moving block technique^ were used to 
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determine the frequency and damping of the response. At least two data points 
were obtained at each test condition and the scatter between the measurements 
was found to be very small. 

Progressing and regressing lag mode data were obtained for the hover test 
conditions, but above about 600rpm the progressing mode was found to be 
contaminated by coupling with the lower stand mode. Because all of the forward 
flight conditions tested were above this rotor speed, no progressing mode data 
were obtained. The regressing mode data should be representative of isolated 
blade results over the entire rotor speed range tested and appears to be of very 
good quality. The hover test envelopes for both the 1982 and 1983 tests are 
shown in Fig. 3(a) and the forward flight test envelopes are shown in Fig. 3(b). 
The edges of the forward flight envelopes were set by the maximum allowable 
rotor flapping loads. 
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Table 1: Measured Model Properties 


Number of Blades 

3 

Radius 

0.81 m 

Chord 

0.0419 m 

Airfoil Section 

NACA 23012 

Non-dimensional Hinge Offset 

0.111 

Blade Inertia About Hinge 

0.01695 kg-m 2 

Blade Mass Center Distance from Hinge 

0 . 188m 

Blade Mass ( Outboard of Hinge) 

0.204 leg 

Non-rotating Flap Frequency 

3.09 Hz 

Non-rotating Lead-lag Frequency 

7.02 Hz 

Average Lead-Lag Structural Damping Ratio 

0.19% critical 

Lock Number y (based on a=5.73) 

7.54 

Lift curve slope a 

5.7 3 

Profile drag coefficient (assumed) 

0.0079 
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5. CORRELATION OF THEORY AND DATA 


We, now come to presenting the correlation between the measured and 
predicted data in hover and forward flight ( 0 H i. 0.55). If not stated 
otherwise, the predicted values based on a linear theory include dynamic inflow 
and uniform steady inflow and the measured values include both the 1982 and. 1983 
test points (figure 3). To facilitate an appreciation of the scope of the 
correlation, we begin with an overview of the data-base. As seen from figure 3, 
the hovering and forward flight cases refer to zero and full structural coupling 
( R=0 and l ) in combination with varying rotational speeds in hover (CKIKIOOO rpm) 
and with two rotational speeds (fl=750 and 1000 rpms ) in forward flight. Further, 
we have four values of collectives in hover( 9 0 = 0”, 4°, 6° and 8“) and three 
values of collectives in forward flight ( 0 O = 0°, 3“ and 6°). Compared with 
hover, the data-base in forward flight is much broader in scope, since for each 
collective, the shaft tilt varies from 0° to 20“ in combination with different 
advance ratios as shown in figure 3b. Thus, given the extensive scope of the 
data-base and the corresponding predicted values, only a small portion of the 
total correlation is presented for the damping levels. (Frequency-data 
correlation is available in reference 17.) We also emphasize how far we can 
isolate different ingredients that participate in the correlation. In 
particular, those ingredients refer to nonuniform steady inflow and unsteady or 
dynamic inflow. In Several cases, we also refer to nonlinear airfoil drag in 
substall ( | a | <.1 2 “ ) and to stall (nonlinear airfoil drag and lift, |a| >12“ ) , and 
heuristicallv identify the deviation between the theory and data with such 
non-linear airfoil characteristics. However, it is good to reiterate the 
following. The nonlinear airfoil section properties including those of drag in 
substall are not included. Though a linear theory is used, we include 
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correlation even when an appreciable portion of the rotor disk (_> 10%) is in 
stall. That inclusion together with stall region contours provides a better 
understanding of the adequacy of the linear theory and also of the possible role 
of airfoil aerodynamic nonlinearity. 

With this as a background, we take up figure 4 which shows the lag re- 
gressing mode damping versus the rotor rotational speed for four values of 
collectives, 0 = 0°, 4“ , 6° and 8 “ . For zero pitch setting, the angle of 
attack is very low. This fact is well reflected in the excellent correlation 
between theory and data except for minor discrepancies for very low rotational 
speeds ( sav, ft < 300 rpm ) . Those discrepancies are consistently observed for 
all the four cases of collectives for ft< 300 rpm, when the predicted damping 
level is less than the data. The reasons for this are not known. Also for 0 O = 
4°, 6" and 8“, the correlations is generally good. However, for ft > 750 rpm, the 
theory deviates from the data and that deviation increases with increasing blade 
pitch (compare the case for 0 O = 4" with that for 0 O * 8“ ) . The 1.62-m model 
was tested in a 2.13 x 3.05 meter test section and the ratio of model height to 
rotor diameter was only 0.63, a relatively low value. Therefore, additional 
damping values were predicted with the empirical corrections for ground and 

recirculation effects (that is, replacing v by kv, where k<l for ground effect 

and k>l for recirculation). It was found that those deviations between theory 
and data are not associated with ground and recirculation effects. Though our 
approach to capture those effects is empirical, recent experimental evidence 
supports this finding. ^ Thus, in summary, figure 4 shows that the data are in 
general agreement with the theory except for ft < 300 for all 0 O , and for ft> 750 
and 0°°>6" . 
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Coming to the forward flight case, we start with the damping data for the 

** 

zero-collective case with £2 * 1000 rpm and R = 0, figures 5 and 6. Figure 5 
shows the damping data as a function of advance ratio for discrete values of 
shaft tilt a s . Figure 6 is a cross plot of the data in Figure 5, that is, 
damping as a function of a s for discrete values of advance ratio p. Thus, 
figure 5 together with figure 6 provides an improved picture of the correlation 
between theory and data for increasing values of shaft tilt a s and advance ratio 
p. Overall, figures 5 and 6 show fairly good correlation and exhibit two main 
features. First, compared to the correlation for high values of pa s , the 
correlation is much better for low values of pa s . For example, at low advance 
ratios ( p 0.15), the correlation is remarkable even quantitatively, virtually 
for the entire range of shaft tilt <0 _> ct s 20 u ). So is the case for low 
values of shaft tilt (a s <_4'') and for advance ratios as high as 0.55, see figure 
5 for ct s = 4°. In fact, right up to the high advance ratio case of 0.3, the 
correlation is equally remarkable for a s <8°. For p=0.4 and a s <4“ (figure 6), 
the correlation is satisfactory although the theory slightly underpredicts 
damping and the reasons for these minor discrepancies of underprediction are not 
known. That the theory shows remarkable agreement with the data for low values 
of pa s is to be expected. After all, the pitch setting is zero and therefore 
the angle of attack is very low for low values of pa g . Coming to the second 
feature, we observe that with increasing p and a s , the theory consistently 
overpredicts damping and increasingly deviates from the data. For example see 
figure 6 for p=*0 . 3 and a s >10° and for p = 0.4 and a s >6' y . This is also to be 
expected, because the angle of attack is found to increase rapidly with 
increasing values of ptana s . 
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The preceding correlation for zero pitch setting shows that the discrepancy 
b’etween theory and data increase with increasing values of ptana s . To 
facilitate explain this discrepancy, we refer to the stall regions of the rotor 
disk based on equilibrium values and use the percent stall area of the rotor 
disk as an objective measure of the stall effects on the theory. Figures 7a and 
7b show two typical stall region (|a|<12°) contours for p=0.3 and 0.5. We also 
observe in passing that for p=0.2 and a s <^20° , the stall region (not shown) 

CLbOUt iZ%- of the rotor dislc.^ This observation corroborates the fact 
that the theory exhibited good agreement with the data for p=0.1 and 0.2 in 
figure 6. As seen from figures 7a and 7 b vnoretha/? 20% of the rotor disk is in 
stall for p=0 . 3 and a s =16 u and for p=0 . 5 and a s =8". It is also seen that 
slightly more than 10% of the rotor disk is in stall for p=0.3 and Ci s =12° 
(figure 7a), for p=0.4 and a s = 6" (reference 17) and for p=0.5 and a s =4° (figure 
7b). Since the blade has zero pitch setting, the area of the stall region 
increases rapidly with increasing values of p and a s , and at the low test 
Reynolds number («170,000) the effects of airfoil aerodynamic nonlinearity seem 
to become more pronounced. Therefore, with increasing values of p and a s , those 
effects seem to mask the predicted data increasingly and account for the 
increasing deviation between theory and data, and further, the prediction 
becomes increasingly suspect. A similar situation exists for the second set of 
data for 3^ collective, taken up next. 

While figure 8 refers to zero flap-lag coupling with £2=1000 rpm, figure 9 
refers to full flap-lag coupling with £2=750 rpm. Both figures essentially show 
the same feature. That is, the correlation is qualitatively accurate only for 
low values of p and ci s . By comparison, the correlation becomes qualitatively 
inaccurate with increasing values of p and a s , the data showing decreasing 
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damping and the theory predicting the opposite trend of increasing damping. To 
% 

fill in details, we consider figure 8. For p=0.2, even for a s * 20“, the stall 
contour results (not shown) show that not more thanf®% of the rotor disk 
experiences stall. ^ Thus, stall is not an issue for p<0.2. However, it is 
seen that there is consistent difference between the theory and data. For 
example at p=0.05 in figure 8, the theory under predicts damping for all values 
of a s . This is possibly due to nonlinear drag in substall at the low test 
Reynolds number. Figures 7c and 7d also show that an appreciable portion of the 
rotor disk experiences stall. For example, for p=0.3 and a s =14 a , nearly 10% of 
the rotor disk is in stall and the stall region for p=0.4 rapidly increases for 
a<.>8“. Thus in summary, nonlinear air foil section characteristics seem to be 
contributing to the discrepancies between theory and data in figures 8 and 9. 

With figure 10 and 11, we come to the data-base for 0° = 6° and p < 0.2, 
also see figure 3b. The data are limited to low advance ratios and stall is not 
an issue here. For the complete range of shaft tilt (0 < a s “<^ 20), the data are 
available only for p = 0.05. The corresponding correlation is given in figures 
11a ( R=0 ) and lib ( R=1 ) . It is seen that the overall prediction is not satis- 
factory. In fact, for ot s *20“ ( 0<ji<0. 15) the prediction shows that damping trend 
increases with increasing p, a trend that is not supported by the data. This 
feature is quite similar to the one referred to earlier in figures 8 and 9 for 
a s >8“. Given the appreciable deviations between the theory and data for the 
hovering case with comparable pitch setting at Q=1000 (e.g. figure 4 at 0“=6°), 
the deviations in figures 10 and 11 for p=0.05 are not unexpected. Inspite of 
the appreciable underprediction, the data supports the theory in showing that 
the damping level virtually remains constant for the entire range of shaft tilt 
at p=0.05. 
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•> We mentioned earlier that the preceding prediction includes uniform steady 
inflow and dynamic inflow. The final set of correlation as presented in figures 
11, 12 and 13, is used to isolate the effects of nonuniform steady inflow and 
dynamic inflow. In figure 11, the predictions from the linear theory are shown 
by full lines. The predictions from that theory refined to include nonuniform 
steady inflow and the predictions from that theory without dynamic inflow are 
respectively shown by broken lines with dots and by broken lines. It is seen 
that nonuniform steady inflow hardly influences correlation. By comparison, 
dynamic inflow improves correlation throughout, but not significantly. Compared 
to the case with R = 1, when the improvement is at best marginal, the im- 
provement is much better for the case with R = 0. As for dynamic inflow, we 
further explore its effects in figure 12 for the hovering case, and in figure 13 
for the forward flight case at high advance ratios of 0.3 and 0.45. Figure 12 is 
a partial reproduction of figure 4 for 8 0 = 6“ and 8". The improvement due to 
including dynamic inflow is noteworthy. However, for 8>800 rpm, the increasing 
discrepancy between theory and data, particularly with increasing pitch setting, 
remains virtually unaffected, even with inclusion of dynamic inflow. The 
correlations in figures 11 and 12 do demonstrate that it is highly desirable to 
include dynamic inflow in predicting inplane damping. The correlation in figure 
13 based on the 1982 data is extremely interesting. Here we have 8 o =0“, 8=1000 
rpm and R=0 at high advance ratios of 0.3 and 0.45. As seen from figures 7a and 
7b, we should expect appreciable stall effects for 0i s >10° at p=0.30, and for 
Ci s >4“ at p=0 . 4 5 , to render the theory suspect. Two features are worth 
mentioning. First, in substall the improvement due to dynamic inflow is 
negligible for p=0.3. For p=0.45 and a s <4° , when the linear theory is reasonably 
valid, the differences in the predicted values without and with dynamic inflow 
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are too small to isolate dynamic inflow effects. Second, the predicted data 
(without inflow), which also does not account for stall shows ’’better cor- 
relation” for the entire range of shaft tilt. However, stall effects increase 
with increasing values of p and a s and mask the prediction without and with 
dynamic inflow. The prediction in figure 13 needs to be resolved for the 
effects of airfoil aerodynamic nonlinearity for a valid treatment of isolating 
dynamic inflow effects. 


6. CONCLUDING REMARKS 

Thus far, we presented the correlation between theory and data on the lag 
regressing mode damping levels. That correlation leads to the following 
remarks: 

1. In hover, the predictions are in general agreement with the data. However, 
some discrepancies at high rotational speeds and blade pitch settings are 
identified. 

2. In forward flight at 0 o =O° and 3°, the correlation is fairly good when 
not more than about 10% of the rotor disk is in stall. When an appreciable 
portion of the rotor disk is in stall (sav 10& or more), the theory is at 
best qualitatively accurate for S^O'' , in that the theory and data show the 
same trend of increasing damping with increasing a s and p, though there is 
appreciable quantitative discrepancy. However, this qualitative accuracy 
is suspect since the nonlinear airfoil characteristics are expected to mask 
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,, the prediction increasingly with increasing values of p and a s . For 8 0 =3°, 

the theory and data respectively show the opposite trends of increasing and 
decreasing damping levels with increasing a s and p. 

3. In forward flight at 0 O = 6 % the data are limited to low advance ratios (p 

<0.20) and virtually the entire rotor disk is in substall. Still, the 
theory significantly deviates from the data. In fact, for os s =20, the 
prediction shows appreciably increasing damping with increasing p, whereas 
the data shows slowly decreasing damping with increasing p. For p=0.05, 
the data confirm the prediction that the damping level essentially remains 
constant for the entire range of shaft tilt ( 0 ^a s ^ 20 °), and the 
discrepancies between the data and prediction are of the same order of 
magnitude as those in hover for 60 = 6 ” and ft >_ 7.50 rpm. 

4. Effects of nonuniform steady inflow hardly affect the predictions. 

5. The inclusion of dynamic inflow consistently improves correlation 
appreciably in hover and marginally in forward flight. However, at high 
advance ratios, cases are identified for which predictions without and with 
dynamic inflow are masked by the effects of airfoil aerodynamic 
nonlinearity to pass a judgement on dynamic inflow effects. 
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Figure 2. Model installed in the aeromechanics laboratory's 7-by-10-foot wind tunnel. 
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Part 2 - Stall Effects 


ABSTRACT 

Experiments have been performed with a 1.62m diameter hingeless rotor in 
a wind tunnel to investigate flap-lag stability of isolated rotors in forward 
flight. The three-bladed rotor model closely approaches the simple theoretical 
concept of a hingeless rotor as a set of rigid, articulated flap-lag blades with 
offset and spring restrained flap and lag hinges. Lag regressing mode stabi- 
lity data was obtained for advance ratios as high as 0.55 for various com- 
binations of collective pitch and shaft angle. The prediction includes quasi- 
steady stall effects on rotor trim and Roquet stability analyses. Correlation 
between data and prediction is presented and is compared with that of an earlier 
study based on a linear theory without stall effects. While the results with 
stall effects show marked differences from the linear theory results, the stall 
theory still falls short of adequate agreement with the experimental data. 

NOMENCLATURE 


a Linear lift curve slope 

c Blade chord 

c-j (a) Lift coefficient at a 

c<j(a) Drag coefficient at a 

ci 0 Lift coefficient at a 

Cd 0 Drag coefficient at a 

c lao Slope of the lift-coefficient curve at a 
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c dao 

d 

e 

F yk 

F zk 

I 

1 

Pi »P2, • • *PlO 
R 
r 
U 

Up 

Ut 
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u To 

a 

a 

°s 

B 

Y 

C 

AUp 

AUy 

Aa 

Bo 

♦ 

P 

n 

W C 

A 

*1 

P 

♦k 


Slope of the drag-coefficient curve at a 
Drag force per unit length of the blade 
Hinge offset/R 

Force per unit length of the k-th blade In the plane of 
rotation 

Force per unit length of the k-th blade In the plane perpen- 
dicular to the plane of rotation 
Blade moment of Inertia 
Lift force per unit length of blade 
Coefficients as defined In equation (13) 

Rotor blade radius or flap-lag structural coupling parameter 

Spanwlse station from hinge 

Resultant flow velocity at a blade sect1on/(QR) 

Transverse velocity component/ (OR) 

Tangential velocity component/ (RR) 

Trim-state transverse velocity component/ (fiR) 

Trim-state tangential velocity component/ (RR) 

Blade section angle of attack 

Trim-state blade section angle of attack 

Shaft-tilt or shaft angle 

Flapping angle of the blade 

Lock number, pacR 4 /I 

Lead-lag angle of the blade 

Perturbed transverse velocity component/ (RR) 

Perturbed tangential velocity component/ (OR) 

Incremental angle of attack, a-a . 

Collective blade pitch 
Local Inflow angle 
Tan-l(U Po /U To) 

Up 0 /U T o 
Air density 
Rotor speed 

Dimensionless rotating lag frequency/R 
Total Inflow ratio, A^ + y tana s 
Induced flow due to thrust 
Advanced ratio 

Azimuth angle of the k-th blade 
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INTRODUCTION 


(’) — < > 
d\p 


In an earlier study* we presented a correlation between data and prediction 
for flap-lag stability In forward flight. The correlation was based on a 
comprehensive data base on the lag regressing mode damping of an Isolated rotor 
with three blades. Virtually rigid flap-lag blades were used, and the data base 
included aerodynamical ly demanding test cases with advance ratio, y, as high as 
0.55 and shaft-tilt angle, a s> as high as 20°. The prediction was based on a 
linear quasi-steady aerodynamics theory with dynamic Inflow. Overall, the pre- 
diction was found to be adequate for very low values of collective pitch, but 
deteriorated at the higher pitch angles. References 2 and 3, though restricted 
to hovering conditions, showed that Inclusion of quasi-steady, nonlinear airfoil 
characteristics (stall effects) significantly Improves correlation. The purpose 
of this continuing study Is to Investigate the effects of nonlinear airfoil charac- 
teristics under more demanding forward flight conditions. 

In forward flight the number of correlation studies based on models that 
are intentionally simplified to Isolate one aspect of the overall rotor 
stability problem are limited. Specifically stated, the structural simplicity 
of a rigid blade model facilitates isolation of aerodynamic effects. For an 
improved picture of airfoil effects, we now investigate the effects of nonlinear 
local lift, of nonlinear local drag and of profile drag at zero angle of 
attack. Furthermore, the correlation Includes different combinations of pitch 
settings (o s 0 O s 8°), advance ratios (o £ y * o*55) and shaft-tilt angles (o 
^ a $ ^ 20°) and thus provides a range of rotor loading conditions In forward 
fl ight. 


EXPERIMENT AND DATA 

For completeness, we Include a brief account of the experimental model, 
for details see reference 1. To ensure the validity of using a simple flap-lag 
analysis for correlations, the three-bladed rotor used flexures to simulate 
articulated blades with spring restraint and coincident flap and lag hinges. 
The effective hinge offset was 0.11R. The blades were stiff relative to the 
flap and lag flexures so that the first flap and lag modes essentially Involve 
only rigid body blade motions. Further, the design Insures a rotating first 
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torsional frequency of at least 9/Rev over the entire rotor speed range 
tested, and thus virtually eliminating the need to consider a torsional degree 
of freedom. The measured and assumed parameter values of the test model are 
given in table 1. 

The model tested was a 1.62m diameter three-bladed hingeless rotor mounted 
on a very stiff rotor stand so that the stability data was representative of an 
isolated rotor, see Fig. la. The rotor had no cyclic-pitch control and collec- 
tive pitch was set manually prior to each run. At a given advance ratio, rotor 
speed and collective pitch the shaft tilt was the only means of controlling the 
rotor. Thus the rotor was operated untrimmed with an unrestricted tilt of the 
tip-path plane. The model was shaken in roll to excite the lag modes, and was 
then locked up and the transient response recorded. The frequency and damping 
data were then obtained from the time histories via the moving block technique. 

In forward flight two rotor speeds were used (Q * 750 and 1000 rpm) corre- 
sponding to values of 0.62 and 0.72. Advance ratio, shaft-tilt angle and 
collective pitch were varied to cover the test envelope shown In figure lb. 

At each condition tested at least two separate damping measurements were obtained. 


Table 1: Model Properties 


Number of blades 

3 

Radius 

0.81 m 

Chord 

0.0419 m 

Airfoil section 

NACA 23012 

Lift curve slope a 

5.73 

Profile drag coefficient at zero 
angle of attack (assumed) 

0.0079 and 0.012 

Nondlmenslonal hinge offset 

0.111 

Blade Inertia about hinge 

0.01695 kg-m 2 

Blade mass center distance from hinge 

0.188m 

Blade mass (Outboard of hinge) 

0.204 kg 

Nonrotating flap frequency 

3.09 Hz 

Nonrotating lead-lag frequency 

7.02 Hz 

Average lead-lag structural damping ratio 

0.185% critical 

Lock number Y (based on a=5 . 73) 

7.54 
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ANALYSIS 


The analytical model consists of an articulated rigid-blade flap-lag model 
with flap and lag spring restraints at the offset hinge. The hinge offset of 
11.1% in flap and lag Is accounted for in both the rotor trim and stability 
analyses. The rotor used no cyclic pitch so that at a given advance ratio, the 
collective pitch and the shaft tilt angles are the known trim parameters. The 
analytical model has the capability to Include stall effects In both the rotor 
trim (zero cyclic pitch, and cyclic flapping present) and Floquet stability 
analyses. At an azimuth angle we consider a blade element of length dr of 
the k-th blade. The corresponding resultant flow velocity U and Inflow angle <|> 
in terms of normal and tangential velocity components are given by 


U = / U T 2 + Up 2 S U T (H4SUp 2 /U T 2 ) 

(1) 

(j> = tan-lUp/Uy 

(2) 


The force components parallel and perpendicular to the plane of rotation are as 
fol lows: 


dFyk * (-1 sin - d cos $)dr 
dFzk = (1 cos $ - d sin $)dr 
where 

1 = HpcQ 2 R z u2 Cl (a) 
d = ttpcfl 2 R 2 U 2 c d (a) 


For a pitch setting 8 0 and trim angle of attack a, we have 


a = 0 O - <j» and a * 8 0 - $ (6) 

where 


♦ = tan-1 Up 0 /U To (7) 


The perturbation angle of attack Aa about a Is expressed as 


a - a = Aa = - <t> (8) 

~ * 11)03 Uf> x Up3 

~ *° " T ' ^ + 3U? 


3(a). 

3(b) 

4(a) 

4(b) 
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Expanding q(a) and Cd(a) about a we get 

C 1 (a) = C] (a) + (^1)- (a - a) (9a) 

c d( Q ) = c d ( a ) + (^)a (a - a) (9b) 

dot ' 

Similarly substituting equations 1, 4 and 9 In equation 3, we get 
dFyk = HpcR 2 n 2 U{ - (Cd 0 + Cd a0 Aa) Uj - (cj 0 + ci ao Aa)Up}dr (10a) 

dF zk = ^pcR 2 fi 2 U{ - (Cd 0 + Cdao Aa ) U P + (ci 0 + ci ao Aa) Uj}dr (10b) 

where 

Up = e cos C + cos 8 ( 1 + C)r/R + p sin ( 4 ) + C) (11) 

Uj = e sin C sin 8 + r8/R + A cos 8 + p sin 8 cos (t + C) 


Now we perturb the total velocity component about the trim state. That is 
Up = U Po + AU T ( 12 ) 

U T ■ u To + AUj 

Substituting equations 1, 8 and 12 In equation 10 and neglecting the products 
of perturbation quantities such as AUj 2 , we get 

dF yk * l$pcQ 2 R 2 U{ pi Uj 2 + P2 U T U P + P3 U P 2 + P4 U P AUp + P 5 Up AUj}dr,(12a) 
dF zk = ^pcfl 2 R 2 U{pg Uj 2 + p 7 UjUp + pg Up 2 + pg Up AUp + pig Up AUj}dr,(12b) 
where 

A 3 

Pi = { - Cdo - Cdaot^o “ ^ — ) 1 (13) 

P2 = I c dao - c l0 - ci ao (^ 0 - -j-) ? 

P 3 = { - 1/2 Cdo - ^ 00 - 1/2 <J> 0 C] 0 + ci ao (l — —) } 

2 

P4 = f c dao (1/6 * 0 - 1/2 * 0 3) - 1/2 * o c l0 - c la0 Ig- ) 
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PS = { - + C '0 ^ + * C1ao+ ° 3 } 

P6 = { c lo + c lao($o " ^ 4>o 3 ) } 

P7 = f " c lao - c do _ c dao ($o ““^o^) i 

p 8 = { £1° + 5 1q °! ° - - C -dolo + c da0 (1 - “ 4> 0 3 ) } 

2 3 2 J 

P9 - ( - c la0 (l/6 - - £dao^ } 

3 2 6 

PlO = { c 1qo^o 2 - c do^o 2 + Cdaol’o 3 } 

6 2 6 

We observe that the quantities q a0 , Cdao, c lo and c do 1n Pi* P2»**»PlO 
experience large variations In numerical values and that the derivation Is 
carried out In terms of Pi,..,Pio throughout (that Is without breaking the p 
quantities In terms of Individual components). Further we stipulate that 
P 1 ....P 10 are of order unity In deriving the equations. The subsequent derlva 
tlon of the flap-lag-dynamic Inflow equation follows reference 1 except for the 
following difference. At any azimuth station, we Integrate the aerodynamic 
terms numerically along blade span since c^q* C]qo, Cdo» Cdao are complex func 
tions of radial coordinate and azimuth. 

To generate the equations, we use a special purpose symbolic processor DEHIM 
(Dynamic Equations of Helicopter Interpretive Models). 6 * 7 It also generates 
FORTRAN coded statements of the equations which are utilized to form subroutines. 
These subroutines are directly linked with numerical analysis program to facili- 
tate evaluate the coefficients of the governing equations. The numerical analy- 
sis program evaluates the rotor trim parameters with stall characteristics, 
performs Floquet stability analysis and Identifies the modes. The numerical 
integration of aerodynamic terms In the equations along the span Is done by a 10 
point Gaussian quadrature. We generate the Floquet transition matrix by 
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subroutine DVERK(IMSL) which Is a Runge-Kutta-Verner fifth and sixth order 
method. The program evaluates the four aerofoil characteristics ci 0 , ci ao , Cd 0 
and ddao f° r an y given angle of attack a by linear interpolation from a table 
with data at 5°1nterva1s for o s a s 360°. The airfoil characteristics used In 

this paper are shown in figure 2. For low angles of attack these are the same 
as the airfoil properties used In our earlier study 1 based on the linear theory 
except for the following difference. In reference 1, c^ (a « 0) » 0.0079 
throughout. However, In the present study we consider two values of cj (a « 0); 
0.0079 as in references 2 and 3 and 0.012 as in references 4. With double pre- 
cision arlthmatlc, the average CPU time for each case was about 5 minutes on VAX 
750 computer. 

Finally, we conclude this section with a note concerning the comparison bet- 
ween the numerical results from the stall theory and from the linear theory. 

For a consistent comparison for the entire test envelope, we compute the damping 
data from the stall theory and also from the same stall theory by suppressing 
the stall effects. The latter computations refer to the linear theory and com- 
pare with those In reference 1. 
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FIG. 2 ASSUMED NON LINEAR LIFT AND DRAG CHARACTERISTICS 
OF AIRFOIL 
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RESULTS 


We predict the lag regressing mode damping for the parameter values of the 
test model given In Table 1, covering the entire test envelope shown In figure 
lb. We use two theories, the quasi-steady aerodynamics with dynamic Inflow 
(linear theory) and the stall theory which Is this linear theory refined to 
Include nonlinear airfoil section local lift and drag characteristics In a 
quasi-steady manner. If not stated otherwise, the following convention and para- 
meter values apply: full lines for predictions with linear theory, dotted lines 
for predictions with stall theory and c^ (a * o) * 0.0079, R * 0.0 and * 
0.72 (ft = 1000 rpm). 

We begin the discussion of correlation with figure 3 which is for the 
hovering case for four values of the collective pitch setting 0 O ° a o, 4, 6 and 
8. The Improvement with Inclusion of quasi-steady stall is marginal and It Is 
due to nonlinear airfoil-section local drag coefficient (non-linear drag for 
short) In substall. This Is expected because the angle of attack Is low for all 
the four cases tested. For example at 0.7R, the approximate mean or trim angle 
of attack a varies from 0° for 0 O ° a 0 , to 4° for 0 O ° = 8. Overall the correla- 
tion is fair. However we observe two types of consistent underpredictions. For 
ft < 300 rpm, the first type is observed for which the deviation from the data 
essentially remains the same with Increasing pitch setting. For ft > 700 rpm, 
the second type Is observed for which the deviation Increases with Increasing 
0 Q . These deviations were found not to be associated with ground effect, recir- 
culation and nonuniform steady Inflow. 1 The present stall theory shows that the 
deviations are not associated with nonlinear drag. The deviations are surpris- 
ing and merit further study. (The role of dynamic Inflow with stall and of 
higher airfoil profile drag with cj (a= 0) = 0.012 Is discussed later). 

We now present the forward flight case In figure 4 which shows the correla- 
tion for zero collective and for relatively low values of the shaft-tilt angle 
(a s °£ 6). For p s 0.4, both the theories show good agreement with the data for 
all the four cases, a s p * 0, 2, 4, and 6. This Is expected since the predic- 
tions refer to low thrust conditions due to zero pitch setting and low shaft- 
tilt angles. For p > 0.4 and a s ° * 4 and 6, the non-linear effects begin to 
affect the predictions. 

To facilitate further discussion, we refer to the areas of the stall re- 
gions (| a | > 12°) based on trim values as a means of quantifying stall effects. 
Stall plots are given In figure 5 Including different combinations of 0 O anc j Qs 
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for various y values. We observe (not Included In figure 5 but given reference 
1) that 10% to 12% of the rotor disk Is in stall for y » .5 and a s ° = 4, and for 
y = .4 and a s ° * 6, and that this percentage Increases to 16-18 for y = .5 and 
a s ° = 6. Overall, both the theories give good correlation. However, a quali- 
tative aspect of the correlation merits special mention. The data point at y * 
0.55 and a s ° = 6 shows that the damping is Increasing with increasing y; a 
trend that Is not captured by the stall theory. We should also mention that 
this trend Is captured by the linear theory, although nearly 1/4 of the rotor 
disk Is experiencing stall for y = .55 and a s ° * 6. 1 

Figure 6 shows the correlation at the same zero collective pitch setting 

for higher values of the shaft-tilt angle (8 s a s ° £ 20) from hovering to high 

advance ratios (o £ y £ .55). Both the theories show close agreement with the 

data for all four cases (a.° * 8, 12, 16, and 20) when less than 10%-12% of the 

1 5 

rotor disk Is In stall. Specifically, these ranges are: y £ .3, a s ° - 8; y 

£ .25, a s ° = 12; y s .225, a s ° = 16, and y £ .175, a s ° ■ 20. They represent 
very low thrust conditions with negligible Influence of nonlinear drag in 
substall. Moreover, the stall theory falls to capture the trend of the data 
( 1 . e . , Increasing damping with increasing y ) when more than 20-25% of the disk 
Is in stall, as observed for the following ranges: y >.425, a s ° ■ 8; y >.325, 

a s ° = 12; y > .3, a s ° * 16, and y = 0.25, a s ° * 20. For the In-between cases 

when 12%-20% of the disk Is In stall, the stall theory shows perhaps a slight 
Improvement. Reiterating we summarize the correlations In figures 4 and 6 as 
follows. First, when more than about 20-25% of the disk Is In stall, which 

occurs for high values y and a s , the stall theory does not capture the trend of 

the data. We suggest that dynamic stall may be contributing to this situ- 

ation. It is not known why the linear theory gives better results here. 
Second, when less than about 10-12% of the disk is in stall, the stall theory Is 
found to merge with the linear theory (low a values due to low values of pitch 
setting, y and a s ). Third, between these two ranges the nonlinear theory seem 
to give slightly better results than the linear theory does, but this Is not 
conclusively demonstrated by the nonlinear theory. 

For the three-degree collective case, figure 7 shows the results for 0 
£ y ^ .55 and 0 s a s ° £ 20. The data shows the trend that the damping decreases 
smoothly with Increasing values of y and a s (a s £ 12°). The linear theory beco- 
mes qualitatively Inaccurate with Increasing y and a s ° (£ 12), showing the oppo- 
site trend of Increasing damping with Increasing y. The stall theory, like the 
data shows an eventual damping reduction with Increasing advance ratio, but the 
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character of this reduction Is quite different. For small values of shaft-tilt 
angle a s ° (s 8), the data Is available for y £ .25 and stall Is hardly an Issue. 
The minor differences between the two theories are due to nonlinear drag In 
substall. That the stall theory shows sudden decrease In damping at y « .4 for 
a s ° = 0 Is likely due to sudden Increase in stall effects, as seen from the 
stall plots In figure 5 for e 0 0 - 3 and a s ° ■ 0. Further, as seen from figure 7 
for a s ° = 0, 4, and 8, the difference between theory and data remains nearly the 
same as It was In hover. According to an earlier study, that difference was 
found not to be associated with nonuniform steady Inflow (more discussion on 
this difference later In figure 13). For higher shaft-tilt angles (a s ° z. 12) 
in figure 7 the stall effects are negligible for the following three ranges of 
data: y £ .275, a s ° « 12; y £.225, a s ° = 16; and y £ .175, and a s ° * 20. That 
Is, less than 10-12% of the disk Is In stall and there Is negligible difference 
between the two theories. For the remaining three ranges of data with a s ° * 
12, stall effects become Increasingly Important with Increasing y and a s . For 
example, at y = .35 and a s ° - 16, nearly 25% of the disk Is In stall. The stall 
theory substantially differs from the linear theory but still does not correlate 
well with the data. 

Figure 8 shows the correlation for the six-degree collective case. Only a 
limited amount of data are available (y £ .15) and stall Is not an Issue here 
(less than 10% of the disk experiences stall for the entire set of data, also 
see figure 5 for B 0 ° « 6 and a s ° = 16). The behaviour of the theories relative 
to each other remains essentially as It was for the 3° collective case. 

The predictions from the stall theory, as presented earlier, Include the 
total stall effects of lift and drag. It Is Instructive to estimate how much 
of that total Is due to nonlinear lift and how much, due to nonlinear drag. 
This question Is addressed In figure 9 for hovering and In figure 10 for forward 
flight. The hovering case has eight degree collective treated earlier In figure 
3, and substall conditions are present throughout (a « 4° at .7R). We have 
shown predictions for four cases — linear lift and drag In combination with 
nonlinear local lift and drag -- as Identified In the figure. Given the 
substall conditions, the nonlinear lift characteristics have no Impact on the 
predictions. As expected, the predictions from the nonlinear 1 Ift-and-drag 
theory (stall theory) merges with those from the linear-lift and nonlinear-drag 
theory. The same Is true of the other two predictions from the nonlinear-lift 
and constant-drag theory and the linear-lift and constant-drag theory (linear 
theory). Thus the relevant Ingredient below stall Is nonlinear drag, 
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although for the specific case In figure 9, Its Impact Is not appreciable. 
Nevertheless, the simplicity of the hovering case lays the ground work for the 
more demanding forward-flight case taken up next. 

Figure 10 shows the predictions for the same four combinations of lift and 
drag characteristics treated for the hovering case and these combinations are 
Identified In figure 10. We have included all the three cases with 0 O ° * 0, 3, 
and 6 as treated earlier (figures 4, 6, 7 and 8). In substall, the predictions 
are as expected. However, when stall becomes an Issue, that Is, when more than 
10-12% of the disk Is In stall (6 0 ° * 0, \i * 0.225, 8 0 ° * 3, p It .225 and 9 0 ° * 
6, p * .3) the trends of the predictions are extremely Interesting, particularly 
of the predictions from the nonlinear-lift and constant-drag theory. It Is seen 
that the key Ingredient Is nonlinear lift (and not nonlinear drag) that affects 
the prediction qualitatively (Increasing or decreasing damping with Increasing 
V)* 

In figure 11, we address the question of how much better Is the stall 
theory with dynamic Inflow when compared to the stall theory without dynamic 
Inflow. In hover dynamic Inflow Improves the correlation consistently and 
throughout. In forward flight the Impact of dynamic Inflow Is negligible and a 
typical example Is shown In figure 11 for 0 O ° = 6 and a s ° * 16. 

Finally, in figures 12 and 13, we discuss the sensitivity of the predic- 
tions to the assumed values of c^ (a = 0) or Cd.mln values. While figure 12 
refers to the hovering conditions for four values of the collective pitch 
setting (as In figure 3), figure 13 refers to the forward flight conditions for 
a typical cross section of the data In forward flight, as In figures 6,7, and 8. 
In substall, the higher value .012 gives better overall correlation both In 
hover and forward flight. This covers all the hovering data and part of the 
data in forward flight as discussed earlier. When more than 10-12% of the disk 
is in stall, we know from figure 11 that the predictions are qualitatively 
affected by nonlinear local lift and not by nonlinear local drag. This Is well 
borne out by the results In figure 13 which shows that the qualitative aspects 
of the prediction In stall are not sensitive to reasonable changes In the 
assumed Cd >m -j n values. 


CONCLUDING REMARKS 

1. When less than 10-12 percent of the disk is stalled the linear and non 
linear theories give nearly the same results, with the nonlinear theory margi- 
nally better. When 10-20% of the disk Is In stall, the nonlinear theory may 



slightly improve the correlation. However only a limited number of data are 
involved in this range and this improvement is not conclusively demonstrated. 

2. For highly stalled cases in forward flight, when more than about 25 percent 
of the disk Is stalled, the stall theory differs markedly from the linear 
theory, but neither theory does well. However, the linear theory seems to be 
"qualitatively accurate" for 8 = 0°. That the linear theory shows the qualita- 
tive accuracy Is peculiar since appreciable stall effects should be expected. 
That the highly stalled, high advance ratio cases are not well predicted by the 
nonlinear theory perhaps Indicates dynamic stall effects. 

3. In hover, the data is in substall, and the difference between the stall 
theory and data increases with Increasing n and 0 O . This difference is not 
associated with nonlinear local drag characteristics. The correlations were 
conducted for >m j n = 0.0079 and 0.012. The latter higher value gives better 
correlation In hover, and In forward flight In substall. Under stall conditions 
of forward flight, the qualitative aspect of the correlation Is not sensitive to 
perturbations in values. 

4. For stall conditions of forward flight, the key Ingredient Is nonlinear 
local lift coefficient which qualitatively changes the prediction when compared 
with the predictions from the linear theory, and from the theory with linear- 
lift and nonlinear-drag characteristics. Further, while the theory with 
nonlinear-lift and constant-drag characteristics Is close to the stall theory, 
the theory with linear-lift and nonl Inear-drag characteristics Is close to the 
linear theory. 

5. Quasi-steady stall theory with dynamic Inflow Improves the correlation 
somewhat In hover but the forward flight results, while qualitatively very dif- 
ferent than the linear theory results, do not overall show Improved correlation 
with the data. 
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